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Extension of the Thomas-Fermi approximation for trapped Bose-Einstein condensates with an 

arbitrary number of atoms 
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By incorporating the zero-point energy contribution we derive simple and accurate extensions of the 
usual Thomas-Fermi (TF) expressions for the ground-state properties of trapped Bose-Einstein conden- 
sates that remain valid for an arbitrary number of atoms in the mean-field regime. Specifically, we obtain 
approximate analytical expressions for the ground-state properties of spherical, cigar-shaped, and disk- 
shaped condensates that reduce to the correct analytical formulas in both the TF and the perturbative 
regimes, and remain valid and accurate in between these two limiting cases. Mean-field quasi- ID and -2D 
condensates appear as simple particular cases of our formulation. The validity of our results is corrobo- 
rated by an independent numerical computation based on the 3D Gross-Pitaevskii equation. 
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The experimental realization of Bose-Einstein conden- 
sates (BECs) of dilute atomic gases confined in optical and 
magnetic traps fTL 0, Si has stimulated great activity in the 
characterization of these quantum systems. Of particular 
interest are the ground-state properties of trapped BECs 
with repulsive interatomic interactions [4]. These proper- 
ties derive from the condensate wave function ip(r) which, 
in the zero-temperature limit, satisfies the stationary Gross- 
Pitaevskii equation (GPE) ^ 



^V 2 + V(r) + gN\ 
2m 



ip = pii(), (1) 



where N is the number of atoms, g = Anffa/m is 
the interaction strength, a is the s-wave scattering length, 
V(r) = ^rn{u)\r 2 _ + u) 2 z 2 ) is the harmonic potential of 
the confining trap, and pi is the chemical potential. 

Only in two limiting cases can Eq. ([T} be solved analyti- 
cally: in the Thomas-Fermi (TF) and perturbative regimes. 
When N is sufficiently large that p S> Hljj_, Huj z , one 
enters the TF regime. In this case the kinetic energy can 
be neglected in comparison with the interaction energy and 
the GPE reduces to a simple algebraic equation. Useful 
analytical expressions can then be obtained for the con- 
densate ground-state properties [4]. In the simple case 
of a spherical trap characterized by an oscillator length 
a r = ^Jh/muj, Eq. £Q leads in the TF limit to 
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moor + gN \ip{r)\ = p, 



< r < R (2) 



where the condensate radius R = w2fjjhu> a r is deter- 
mined from the condition |^(r)| 2 > 0, and the chemical 

potential p = | (15Na/a r ) 2 ^ 5 fmo follows from the nor- 
malization ofip(r). 

In the opposite limit, when A?" is small enough that the 
interaction energy can be treated as a weak perturbation, 
one enters the (ideal gas) perturbative regime. In this case, 
to the lowest order, ip(r) is given by the harmonic oscillator 



ground state, ip(r) = {na 2 ) 
chemical potential satisfies 



-3/4 



exp(— r 2 /2a 2 ), and the 



(3/2)faj + gn = p 



(3) 



where n = N/ (v27ra r ) 3 is the mean atom density. Away 
from these two limiting cases, in principle, one has to solve 
the GPE numerically. Very few theoretical works have ad- 
dressed the question of looking for approximate analyti- 
cal solutions valid in between the two analytically solvable 
regimes. The most relevant proposals are based on a varia- 
tional trial wave function [6], or on the semiclassical limit 
of the Wigner phase-space distribution function of the con- 
densate [7]. However, the practical usefulness of these ap- 
proaches turns out to be somewhat limited in comparison 
with the simple TF approximation. 

In this work we address the above question from a dif- 
ferent point of view. We start from the usual TF approxi- 
mation and modify it conveniently to account, in a simple 
manner, for the zero-point energy contribution. This en- 
ables us to derive simple and accurate extensions of the 
TF expressions that remain valid for an arbitrary number 
of atoms in the mean-field regime. Specifically, we obtain 
general analytical expressions for the ground-state prop- 
erties of spherical, cigar-shaped, and disk-shaped conden- 
sates that reduce to the correct analytical formulas in both 
the TF and the perturbative regimes, and remain valid and 
accurate in between these two limiting cases. 

We begin by considering a BEC in a spherical trap. In 
principle, we start from the TF relation of Eq. (0. How- 
ever, since we intend to apply this equation to arbitrarily 
small condensates, we introduce a lower cutoff radius r , 
defined through \rmo 2 r 2 = ^Hlj, in order to be consis- 
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tent with the fact that the contribution from the harmonic 
oscillator energy cannot be smaller than the zero-point en- 
ergy. As for the small volume V ~ a 3 corresponding to 
f < r , we do not aspire to get a precise knowledge of the 
wave function therein. Instead, we content ourselves with 
an effective condensate density n Q in that region. As we 
shall see, this is all that is needed to obtain very accurate 
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approximate formulas for most of the condensate ground- 
state properties. Thus we start from the ansatz 
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-mu r 



gN\1>(r)\' 



r < r < R (4a) 



-frw + gy/6/7rn 



fi, < r < r (4b) 



with ip{r) = for r > R. A renormalization constant 
= v/6/7r has been introduced in Eq. ( |4bl to guaran- 
tee the correct perturbative limit. In this limit fx — * \hw 

and R —> r = y3a r . Under these circumstances, only 
Eq. d4bl contributes significantly to the chemical poten- 
tial, and in this case n = N/V . This corresponds to a 
uniform spherical condensate, defined in the finite volume 
V . In order for this uniform density to produce the same 
chemical potential as the ground state of the harmonic os- 
cillator over the volume of the entire space it is only neces- 
sary to renormalize the corresponding interaction strength 
by multiplying by \ZQ/tt. Equations I© also yield the cor- 
rect result in the TF regime. This is mainly a consequence 
of the direct relation existing between the number of parti- 
cles and the size of a trapped BEC. For large condensates, 
such that fx 3> hw, one has R S> r and, as a result, the 
relative contribution from Eq. (l4bl to the normalization in- 
tegral that determines fx becomes negligible. Since we have 
renounced an explicit expression for ip(r) in V , in this re- 
spect, our approach cannot provide more information than 
the TF approach. Only when R^> r can we have a suffi- 
ciently precise knowledge of the wave function and, in this 
case, it coincides with the TF wave function. 

The chemical potential follows from the normalization 
of ip{r). After a straightforward calculation one obtains 

1—5 v3 . .—2 3\/3 / 3\ „ T a 
_* +_(«-!)* -_ ^--j =J V-,(5) 

— , — 2 

where R = R/a r , and fx = -zR is the chemical poten- 
tial in units of hio. As Eq. © shows, the ground-state 
properties depend on the sole parameter Xo = Na/a r . 
When Xo ^ 1 (TF limit) the above equation leads to 
fx = \ (15xo) 2 ^ 5 > as expected. The Xo <C 1 limit cor- 
responds to the perturbative regime and, in this case, one 
obtains Jx = 3/2 + \/2/ir Xo> which is nothing but the per- 
turbative result d3}. For arbitrary Xo> i n principle one has 
to solve numerically the above quintic polynomial equation 
(which has only one physically meaningful real solution). 
This is a simple task that can be carried out with standard 
mathematical software packages. We have found, however, 
a rather accurate approximate solution. It can be shown that 
the expression 
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FIG. 1 : (Color online) Theoretical prediction for the ground-state 
properties of spherical condensates (solid lines). The open circles 
are the exact numerical results. For comparison purposes we have 
also included the TF prediction (dashed lines). 



satisfies Eq. ([5]) with a residual error ifstl smaller than 0.7% 
for any Xo £ [0, oo). Figures QJa) andQJb) show, respec- 

— i — 2 

tively, the predicted chemical potential, fx =5^ , and con- 
densate radius, obtained from Eq. (O (solid lines), along 
with the exact results obtained from the numerical solu- 
tion of the 3D GPE (open circles). For the numerical cal- 
culation we have defined the radius through the condition 
\ip(R)\ 2 = 0.05 |^(0)| 2 . With this definition, Eq. © re- 
produces the numerical R with a relative error smaller than 
3% for any Xo- Most of the error, however, comes from 
the region where Xo 3> 1 (TF limit) because in that region 
R — > Rtf an d it rather satisfies \ip(R)\ 2 = 0. The accu- 
racy with respect to the numerical fx is better than 0.5%. 

A straightforward calculation yields the mean-field inter- 
action energy per particle, e int = e int /Hu> = E int /Nhw, 
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= sxo" 

-6\/3 l h 
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R 7 + V3 (k - 1) R 4 



9^3 L - I 



(7) 



(2/7) M . 



For Xo 3* 1, one recovers the TF result, e in 
In the Xo < 1 limit, using that R = 3 + 2y / 2/Vxo - 
(l/7r)(l/9+A/2/37r)xo+0(Xo) is a perturbative solution 
of Eq. (0, one obtains ej nt = Xo^-'/v / 2vr = gn/2, which 
again is the correct result. Finally, the kinetic and potential 
energies can be readily obtained in terms of the previous 
results by using the exact relations fH 



£kin = Ekin/N - 

E pot /N 



'pot 



fx/2 - (7/4)E int /N, (8a) 
= fi/2 - (l/4)E int /N. (8b) 



I In Fig. [Ha) we show the theoretical prediction for e int , e^ n , 
I and e pot , obtained from Eqs. ©-([8} (solid lines), along 

(6) with the exact numerical results (open circles). 



Next we consider a BEC confined in a cigar-shaped mag 
netic trap with oscillator lengths a 

a z = \/KJrrujj 

2. We shall restrict ourselves to the mean-field regime 



^h/mu)± and 
and an aspect ratio A = uj z /uj ± <C 



which requires NXa\/a 2 > 1 H QU [H. As be- 
fore, we start from the usual TF expression, which we as- 
sume to be valid up to a minimum radial distance r° = 
\/2a±, determined from the condition that the contri- 
bution from the radial harmonic oscillator energy should 
not be smaller than Huj±. This defines an outer region 
V+ = {{r ± ,z): r 2 ± /R 2 + z 2 /Z| F < 1 A r x > r°}, 
which is nothing but the usual TF ellipsoidal density cloud, 
truncated at rj_ = r°. Note that unlike what happens 
with the condensate radius R = \j2[i/huj± a±, which 
remains the same, now the axial condensate half-length 

coincides with the TF 
value Z TF only in the limit n/fux)^ S> 1. For large 
condensates, when [i ^> hw± (TF regime), this is the 
only region that contributes significantly. On the con- 
trary, in the perturbative regime, as /i —* fkj^ most of 
the contribution comes from the inner cylinder V_ = 
{(rj_,z): r± < r° A \z\ < Z}. In this case, if a <C 
aj_, the transverse dynamics becomes frozen in the radial 
ground state of the harmonic trap and the condensate wave 
function can be factorized as ip(r±, z) = (p(r±)<f)(z), with 
l p( r ±) = (vra^) -1 / 2 exp(— r 2 L /2a 2 L ). This corresponds \f\Z 
to a mean-field quasi-lD condensate. Substituting then in 
Eq. £[) and integrating out the radial dynamics, one finds 



l 



2 2 

muj,z 



g 1D N\4>{z)\ 2 = fi, 



(9) 



where g 1D = g/2ira 2 L [12], and we have used that /j, ~ 
huj ± S> hfaw z to neglect the axial kinetic energy. Note 
that g 1D can be conveniently rewritten as gn 2 with n 2 = 
1 / 7r(r° ) 2 , indicating that one can account for the contribu- 
tion from the radial ground state by using a uniform mean 
density per unit area normalized to unity in V_ . Guided by 
these simple ideas, we then propose the following ansatz: 



^muj\r 2 L + ^muj 2 z 2 + gN \ip(r±, z)f = fi, r G V + 
1 



fkjj ± + ^mu} 2 z z 2 + gNn 2 \(p{z)\ 2 



H, r e V- 

with t/j = elsewhere. The normalization of t/j leads to 



^(Vxzf + liVxzy 



NX^, 



(10) 



where Z = Z/a z and R = R/a±. The chemical poten- 
tial 71 = [i/hw± is given by Jl = 1 + ^iy\/XZ) 2 . Now 
the relevant parameter determining the ground-state prop- 
erties is Xi = NXa/a±. When Xi S> 1 (TF regime), Eq. 
CTJJ leads to Jl = ±(15xi) 2/5 and Z = A- 1 / 2 (15xi) 1/5 . 
When *C 1 (mean-field quasi- ID regime), one obtains 




FIG. 2: (Color online) Theoretical prediction for the ground-state 
properties of arbitrary cigar-shaped condensates with A <g 2 
(solid lines). The open circles are exact numerical results ob- 
tained with A = 0.2. The dashed line is the TF prediction. 



/j= l + I(3 Xl )2/3 andZ = A~ 1/2 (3xi) 1/3 , in agreement 
with previous results ifloLfTlll . In general, for arbitrary Xi, 
an approximate solution satisfying Eq. (flOl) with a residual 
error less than 0.75% for any Xi S [0, oo) is given by 



+ 



+ 



(15 Xl )5 + i 57xi + 345 (3*1 



The mean-field interaction energy e int = e int / Hlu ± is 



(11) 
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xzy 



+ l -{Vxz)' 



(12) 



15 X i 

For Xi 3> 1, Eq. (fl~2l) reduces to e int = (2/7)/i, while for 
Xi "C 1, it leads to e int = (2/5)(/x — fouo±), which again 
are the correct analytical limits. As for the condensate den- 
sity per unit length, n 1 (z) = N J 2irr±dr± \ip(r±, z)\ , 
after a straightforward calculation one finds 



ni(z) 



(Vxzy 



1 



z 

~Z 2 



+ 



(Vxzy 



i 



z 2 J 



4a V z ) 16 « 

(13) 

The first term is the contribution from V_ and thus it is the 
only one that contributes significantly in the Xi "C 1 limit. 
On the contrary, the second term, which is the contribution 
from V + , gives the dominant contribution in the Yi S> 1 
limit, in good agreement with previous results ifllTI . 

Figure |2] shows the theoretical predictions for the 
ground-state properties of arbitrary cigar-shaped conden- 
sates with A <C 2, obtained from Eqs. (JU) and <fTTT)— (IT3T) 
(solid lines), along with exact numerical results (open cir- 
cles). 

Finally, we consider a BEC in a disk-shaped trap with 
A > 2 and a z ^> a. In this case, in the mean-field petur- 
bative regime, which occurs when X2 = Na/X 2 a z <C 1, 
the system reduces to a quasi-2D condensate satisfying 
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-mulrl + g 2D N\<p(r ± )\' 



Mi 



(14) 



with g 2 u = g/v2ira z 1 1 3fl - We then rewrite g 2D as 
gn^fii, where fix = l/2a z is a uniform mean density 
per unit length and k 2 x = a/2/vt is the appropriate renor- 
malization factor, and propose the following ansatz: 



^muj 2 z z 2 + ^muj 2 ± r 2 ± + #iV |^(r ± , z)| 2 = /i, 
11 

-hjJ z + -mw^ + gn~ 1 Nn 1 \ip(r ± )\ 2 = /i, 



r£V+ 
(15) 



with t/j = elsewhere. In the above equations, 

{(r ± , z) : ri /i?| F + z 2 /Z 2 < 1 A \z\ > z Q ] and VL 
z): r± < i? A |z| < z }> where z = a 2> -^tf 

Z 



a±, and 



x/2fi/hw z a z . More precisely, one expects k 2 1 — > 
\j2jix in the perturbative regime (^2 <C 1), while k^ 1 — > 
1 in the TF regime (% 2 3> 1). The final results are not very 



sensitive to the specific functional form of k 2 
propose one of the simplest possibilities: 



We thus 



k^( X 2) = V2A + e>(x2-o.i) 

Rtf(X2 



l-y/2fr) 1 



0.1" 



Rtf(X2 

where 0(x) is the Heaviside function and Rtf(X2) 

(15X2)"" 
yields 



(16) 



1/5, 



is the TF radius. The normalization of i\) 
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— Z + -(k 2 
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A 2 ou 



(17) 



where Z = Z/o„ i? = R/a ± , and - R 2 /X = 1. The 

chemical potential is 71 = fJ-/huj z = |(l + i? /A). 

For X2 1 Eq. ( fTTl) leads to the usual TF results, while 
for % 2 <C 1 (mean-field quasi-2D regime), one obtains 7J = 

1/2 + (2 v / l7^x 2 ) 1/2 andii = X 1 ' 2 {&s/2/^ X 2 ) 1/4 ■ An 
approximate solution that satisfies Eq. ( TTtT > with a residual 
error less than 0.95% for any \2 & [0, oo) is given by 

-1/8 



R/V\ = (l/15x 2 ) 8/5 + (%/8x 2 ) 2 



(18) 



After some calculation one finds the following expres- 
sions for the mean-field interaction energy e int = e mt / hw z 
and the condensate density per unit area n 2 (r±): 
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15 
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(19) 
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FIG. 3: (Color online) Theoretical prediction for the ground-state 
properties of arbitrary disk-shaped condensates with A » 2 (solid 
lines). The open circles are exact numerical results obtained with 
A = 20. The dashed line is the TF prediction. 



n 2 (rj 



^(rj-l] , [27J z (r ± )] 3/2 -l 



47raa, 



+ 



6iraa z 



(20) 



where f = (k 2 - 1) and 27J z (r ± ) = 1 + r\{1 - r\/ R 2 ). 

In Fig. [3]we show the ground-state properties of arbitrary 
disk-shaped condensates with A > 2, obtained from our 
analytical formulas [Eqs. dSJ and (fl"8"T l-(l2"0"lll (solid lines), 
along with exact numerical results (open circles). 

In conclusion, modifying the usual TF approximation 
conveniently to account for the zero-point energy contri- 
bution, we have derived general analytical expressions for 
the ground-state properties of spherical, cigar-shaped, and 
disk-shaped condensates that reduce to the correct analyti- 
cal formulas in both the TF and the mean-field perturbative 
regimes and remain valid and accurate in between these 
two limiting cases. Mean-field quasi-lD and -2D conden- 
sates appear as simple particular cases of our formulation. 
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